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Abstract: The effect of vacuum polarization on the propagation of photons in 
curved spacetime is studied in scalar QED. A compact formula is given for the full 
frequency dependence of the refractive index for any background in terms of the 
Van Vleck-Morette matrix for its Penrose limit and it is shown how the superluminal 
propagation found in the low-energy effective action is reconciled with causality. The 
geometry of null geodesic congruences is found to imply a novel analytic structure 
for the refractive index and Green functions of QED in curved spacetime, which 
preserves their causal nature but violates familiar axioms of 5*-matrix theory and 
dispersion relations. The general formalism is illustrated in a number of examples, 
in some of which it is found that the refractive index develops a negative imaginary 
part, implying an amplification of photons as an electromagnetic wave propagates 
through curved spacetime. 



1. Introduction 



Quantum field theory in curved spacetime has proved to be a rich field exhibiting 
many subtle and counter-intuitive phenomena. The most famous, of course, is the 
prediction of Hawking radiation from black holes [1], which has forced a critical anal- 
ysis of unitarity in spacetimes with horizons. More recently, the insights associated 
with holography [2,3] have led to a re-appraisal of the role of locality at a fundamen- 
tal level. Another remarkable, but less well-known, phenomenon discovered during 
the early investigations of QFT in curved spacetime is the apparent superluminal 
propagation of photons due to vacuum polarization in QED. This clearly raises the 
question of whether causality may be violated by quantum effects in curved space- 
time. 

The original result, due to Drummond and Hathrell [4], was obtained by con- 
structing the effective action for QED in a curved background and shows that the 
low-frequency limit of the phase velocity v p h(uj) = c/n(u), where n{uj) is the re- 
fractive index, can exceed the fundamental speed-of-light constant c. This is not 
immediately paradoxical 1 since the 'speed of light' relevant for causality is not v p h(0) 
but the wavefront velocity, which can be identified with the high-frequency limit 
v p h(oo) [12]. In order to settle the question of causality, it is therefore necessary 
to go beyond the low-energy effective action and show explicitly that n(oo) = 1. 
However, a serious problem then arises because of the Kramers- Kronig dispersion 
relation [13-15], which is proved in Minkowski spacetime on the basis of apparently 
fundamental axioms, especially micro-causality, together with standard analyticity 
properties of QFT amplitudes. This states: 



Since unitarity, in the form of the optical theorem, normally implies that Imri(w) 
is positive, a superluminal n(0) < 1 would seem to imply a superluminal wavefront 
velocity, n(oo) < 1, with the associated violation of causality. 

The resolution of this apparent paradox was found in our recent papers [16, 17]. 
We showed there that generic geometrical properties of null geodesies in curved 
spacetime imply a novel analytic structure for the refractive index which invalidates 
the Kramers-Kronig relation, at least in the form ( |1 . 1| ) . The complete frequency 
dependence of the refractive index was found in simple examples and it was shown 
explicitly how a superluminal n(0) is reconciled with n(oo) = 1, ensuring causality. 

^or a review of the issues involved in reconciling superluminal propagation with causality for 
QED in curved spacetime, sec refs. [5 8] . Related work on causality can be found, e.g., in refs. [9-11] . 




(Minkowski) 



(1.1) 
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An important implication of this result is that the conventional assumptions 
about the analytic structure of amplitudes in QFT, which underpin the whole of 
S'-matrix theory and dispersion relations, have to be reassessed in curved spacetime. 
Because of the intimate relation of analyticity and causality, this is key issue both 
for QFT in curved spacetime and, most likely, for quantum gravity itself. It also 
highlights the danger in theories involving gravity of relying on identities and intu- 
ition derived from conventional dispersion relations to extrapolate from low-energy 
effective field theories to their UV completions. In particular, the occurrence of 'su- 
perluminal' behaviour in a low-energy theory does not necessarily mean that such 
theories do not have consistent, causal UV completions, either in QFT or string 
theory [18-21]. 

The centrepiece of the present paper is the derivation of a formula for the full 
frequency dependence of the refractive index for QED in an arbitrary curved space- 
time expressed entirely geometrically, specifically in terms of the Van Vleck-Morette 
(VVM) matrix in the Penrose limit [22-24]. The calculation uses conventional QED 
Feynman diagram methods with the heat kernel/proper time formulation of the prop- 
agators, rather than the worldline method used in our earlier papers. This allows 
us to retain the critical insight of the worldline approach in motivating the impor- 
tance of the Penrose limit, while strengthening the contact with the well-developed 
differential geometry of null geodesic congruences in general relativity. 

This geometry plays a central role here both in the derivation of our key formula 
for the refractive index and in the interpretation of its analytic structure. In particu- 
lar, the idea of the Penrose limit is vital in establishing the generality of our results. 
The fundamental insight provided by the worldline analysis is that to leading order 
in R\ 2 C (where R is a typical curvature and A c = 1/m, the Compton wavelength of 
the electron, sets the quantum scale), the one-loop corrections to photon propagation 
are governed by fluctuations around the null geodesic describing the classical photon 
trajectory. It is precisely this geometry of geodesic deviation in the original curved 
spacetime that is encoded in the Penrose plane- wave limit. This also explains why 
the final result for the refractive index can be expressed purely in terms of the VVM 
matrix since, as we explore here in some detail, this is in turn determined by the 
Jacobi fields characterizing geodesic deviation [25-27]. 

A crucial feature of the geometry of null geodesic congruences is the occurrence 
of conjugate points, i.e. two points on a null geodesic which can be joined by an 
infinitesimal deformation of the original geodesic [25,27]. Their occurrence is generic 
given the validity of the null energy condition, which is an important assumption 
in most theorems involving causality, horizons and singularities in general relativity. 
The existence of conjugate points implies singularities at the corresponding points in 
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the VVM matrix. Translated into the quantum field theory, these imply singularities 
in the refractive index in the complex tu-plane - in particular, n(oS) must be defined on 
a physical sheet with cuts running on the real axis from to ±00. This novel analytic 
structure has important consequences, most notably the loss of the fundamental S- 
matrix property of real analyticity for the refractive index, i.e. n(uj*) = n(uj)*, which 
is assumed in the derivation of the Kramers-Kronig relation ( |1 . 1| ) . The relation with 
causality means that analyticity is a key property of QFT amplitudes, so we must 
show how, despite the violation of the Kramers-Kronig relation, the new analytic 
structure of the refractive index is reconciled with, and indeed essential for, causality 

It is important to emphasize that the essential physics underlying this discussion 
is much more general than the specific application to the refractive index in QED. 
It shows how the geometry of curved spacetime can modify the analytic structure of 
Green functions and scattering amplitudes in quantum field theory in a quite radical 
way. This is sure to have important physical consequences which we have only just 
begun to explore. Certainly, the implications for S'-matrix theory and dispersion 
relations appear to be far-reaching. 

Of course, a full discussion of causality, and micro- causality, must be framed 
more generally in terms of the Green functions of the theory. In this paper, we 
explicitly construct the one-loop corrected Green functions for QED in the Penrose 
plane-wave spacetime, which is sufficient to address the issues of causality in photon 
propagation. The full range of Green functions — Feynman, Wightman, retarded and 
advanced, commutator (Pauli- Jordan or Schwinger) — is found and they are shown to 
exhibit the expected good causality properties. In particular, the retarded (advanced) 
Green functions are shown to have support only on or inside the forward (backward) 
light cone. This confirms that, even at one-loop, the commutator function vanishes 
outside the light cone, which is the conventional quantum field theoretic definition 
of micro-causality. 

The paper is organized as follows. We begin with a review of the classical theory 
of wave propagation in curved spacetime in the eikonal approximation. We then 
summarize the low-energy effective field theory, extending the Drummond-Hathrell 
result to scalar QED. 2 Our main result is presented in Section 4, where we calculate 
the one-loop vacuum polarization in the Penrose plane-wave limit and derive the 
fundamental formula for the refractive index in terms of the VVM matrix. Section 
5 reviews the geometry of geodesic deviation and a number of important identities 
relating the VVM matrix, geodesic interval and Jacobi fields are derived. The Ray- 



2 The case of spinor QED is similar in terms of physics to the results presented here. However, 
the formalism requires further technical developments and will be presented separately. 
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choudhuri equations are used to demonstrate the generic nature of conjugate points. 

The analytic structure of the refractive index is studied in Section 6. The argu- 
ment leading from the existence of singularities in the VVM matrix to the definition 
of the physical sheet for the refractive index in the complex cu-plane is explored and 
the consequences for the Kramers-Kronig relation and causality are carefully dis- 
cussed. The explicit construction of the retarded, advanced and commutator Green 
functions at one-loop, demonstrating that they have the required causal properties, 
is presented in Section 9. 

These formal results are illustrated in Sections 7 and 8 in a number of examples, 
demonstrating explicitly the predicted relation of the geometry with analyticity and 
causality. The cases of conformally flat and Ricci flat symmetric plane waves, and 
also weak gravitational waves, are calculated in detail, with the latter two exhibiting 
gravitational birefringence. Remarkably, we also find that the refractive index may 
develop a negative imaginary part, contrary to the conventional flat-spacetime ex- 
pectation based on unitarity and the optical theorem. Although the physical origin 
of this effect remains to be fully understood, it corresponds to a quantum mechanical 
amplification of the electromagnetic wave as it passes through the curved background 
spacetime, over and above the geometric effects of focusing or defocusing, apparently 
due to the emission of photons induced by the interaction with the background field. 
Finally, our conclusions are summarized in Section 10. 



2. Classical Photon Propagation in Curved Spacetime 



The classical propagation of photons in curved spacetime is governed by the covariant 
Maxwell equation, 

V M F^ = 0, F^ U = V,A U -V U A^. (2.1) 

In a general background spacetime, it is not possible to solve these equations exactly. 
However, we will work in the eikonal, or WKB, approximation which is valid when the 
frequency is much greater than the scale over which the curvature varies, u ^> VR. 
(Here, R is a measure of the curvature scale, for instance a typical element of the 
Riemann tensor.) In this case, we can write the electromagnetic field in the form 

A,(x) = e,(x)e-^ , (2.2) 

where the eikonal phase 6 is 0{oS) and is O(u ). Substituting into Maxwell's 
equation, and expanding in powers of we find the the leading and next-to- 
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leading order terms are 



e~ ie . (2.3) 



- dQ ■ dQe u + 2idQ ■ V e v + ie"V ■ dQ + ■ 

The leading order term yields the eikonal equation 

8Q-dQ = 0; (2.4) 

so the gradient k^ 1 = d^Q is a null vector field. This vector field defines a null 
congruence, that is a family of null geodesies whose tangent vectors are identified 
with the vector field k^. This vector can also be identified with the 4-momentum of 
photons: in this sense the eikonal approximation is the limit of classical ray optics. 

It is convenient to introduce a set of coordinates (u, V, Y a ), a = 1,2, the Rosen 
coordinates, that are specifically adapted to the null congruence: u is the affine 
parameter along the geodesies; V is the associated null coordinate so that 

e = uV; (2.5) 

while Y a are two orthogonal space-like coordinates. As explained in ref . [23] , the full 
metric g^ can always be brought into the form 

ds 2 = -2du dV + C{u, V, Y a )dV 2 + 2C a {u, V, Y b )dY a dV + C ab {u, V, Y c )dY a dY b . 

(2.6) 

The null congruence has a simple description as the set of curves (u, V, Y a ) for fixed 
values of the transverse coordinates (V,Y a ). It should not be surprising that the 
Rosen coordinates are singular at the caustics of the congruence, that is points 
where members of the congruence intersect. 



The next-to-leading order in the eikonal approximation ( |2.3|) gives an equation 
for the evolution of along a null geodesic: 

k- Ve / * = -|e"V- k . (2.7) 

It is useful to make the decomposition = As^, where i^ 1 is the unit normalized 
polarization vector and A represents the amplitude. Eq. Q2.7|) is then equivalent to 
the two equations 

2.8 

k ■ V log .A = — |V • k . 

At this point, we fix the gauge by choosing A u = along with the condition V^A 11 = 
0. The latter implies the transverse condition k ■ e = whilst the former means 
that we set the component of e along k to zero. Hence, there are two independent 



- 5 - 



solutions for the polarization vector irq, i = 1,2, which we normalize as em = 5ij. 
These span the directions associated to the space- like coordinates Y a . The second of 
eqs.( |2.8| ) relates the change of the amplitude along a null geodesic to the expansion 
9 = V^/^, one of the optical scalars appearing in the Raychoudhuri equations (see 
Section 5). 

Later, when we calculate the one-loop correction to the mass-shell condition we 
shall have to take the photon wavefunctions off-shell at tree level. This can be done 
conveniently by modifying the eikonal phase to 

e=cj(V-# ij (u]w)) . (2.9) 

We have indicated the polarization dependence explicitly, in which case the phase 
can be thought of as a 2 x 2 matrix with 

Ajf = gM e -to(v-*«(«^)) ; (2.10) 

with an implicit sum over j = 1,2, in which case, 

V^ (i) - 2u — £ (i) e , (2.11) 

to leading order in the eikonal approximation. If , & i j(u;uS) is perturbatively small 
then the local phase velocity matrix is c(5ij — ddij(u;uj) / du), which gives a matrix 
of refractive indices 

n(u;uj) = 1 + 2 — )-!—!- . 2.12 
ou 

Notice that in order for the correction to remain perturbatively small, the refractive 
index should strictly-speaking approach 1 in the infinite past and future. In other 
words the spacetime should become flat in these limits. 



3. Effective Action and Low-Frequency Propagation 

The low-frequency limit of the phase velocity, which exhibits the superluminal effect, 
can be found by considering the modifications to the Maxwell equation following 
from the leading terms in a derivative expansion of the one-loop effective action. 
This was the approach taken in the original work of Drummond and Hathrell [4] . 

The generalization of the QED effective action to all orders in derivatives was 
subsequently given in ref. [28,29], extracting the relevant U RFF" terms from the 
general heat kernel results of Barvinsky et al. [30]. (See also refs. [31-34] for re- 
lated heat kernel results.) Although these results were given for spinor QED, it is 
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straightforward to find the corresponding results for scalar QED from the formulae 
in [28]. In particular, this allows us to write the leading-order effective action for 
scalar QED and deduce the corresponding low-frequency phase velocity, providing a 
useful consistency check on our general result for the full refractive index in scalar 
QED. 



The relevant terms in the effective action to one loop are 

T = J d A x y/g[ - \F^F^ + aRF^ + bR^F^F" 
+ cR^xpF^F^ + dW^ x W v F\ + ■ ■ 
where, in the notation of [28], 



(3.1) 



87rm 2 

"■ (h 2 (0) - 2h' (0)) 



(3.2) 



87T771 

d=- -2^(0) . 



57rm 2 



For scalar QED JV = 1 while for spinor QED jV = 2 2 = 4. Notice that we have 
used the identity 

J d A x y/gF^UF^ = J d 4 x^[2V,F^ x V u F\ - 2R pv F» x F\ - R pvXp F» v F x ? 

(3.3) 

to write the action in the form 



The all-orders effective action derived in [28] is expressed in terms of RFF-type 
operators acted on by functions of the Laplacian given in terms of the form factors 
/(□) and F(n\i, D 2 , D 3 ) computed in ref. [30]. In particular, the quantities h(0) can 
be expressed in terms of this collection of form factors as follows: 



h(o) = 1*1(0) - £f s (o) , 

h 2 (0) =F 8 (0) , 

/ i3 (0) = -iF 3 (0) + F (0) . 



(3.4) 



For spinor QED [28] 



/Ko) = -£, &{o) = -m> F M = h 

Fs(0) = , F 3 (0) = i , F o (0) = i . 



(3.5) 



-7- 



So the coefficients are 

a 13a a a 

a = , b = , c = , a = . (o.v) 

lAAum 2 360nm 2 3607im z SOrcm 2 

reproducing the original Drummond-Hathrell effective action [4]. 

For scalar QED, we can readily see that f\ = 0, F± = and F% = 0, while the 
other quantities are as above, so the coefficients in this case are 

ot ot ot 

a = , b = , c = , d = . (3.7) 

720vrm 2 UAOirm 2 480vrm 2 v ; 



To relate the effective action to the calculation of the refractive index we write the 
modified Maxwell equation corresponding to (|3.1|) and substitute the eikonal ansatz 
fl2.10|) . This gives the general result for the low-frequency limit of the refractive 
index [4, 7] 

; 0) = dij + 2bR uu (x)5ij + 8cR uiuj (x) . (3.8) 

Notice that the refractive index, and the phase velocity, is a local quantity in space- 
time. Although we call this the "low-frequency" limit, we are still working in the 
eikonal approximation. Low frequency refers to the fact that the dimensionless ratio 
uj\fR/m 2 is small. With this in mind, for spinor QED we find 

riij(x; 0) = Sij - ^ 2 (l3R uu {x)8i j + AR uiuj (x)) , (3.9) 

while for scalar QED, 

n i:j (x; 0) = Sij + 3g Q " m 2 { R nn{x)S ij - 2R uiU j(x)) . (3.10) 

Notice that the opposite sign of the h coefficient means that scalars and spinors 
respond oppositely to the Ricci curvature. Since the null energy condition requires 
Ruu(x) > 0, this means the low-frequency phase velocity is superluminal for spinors 
in a conformally flat background, but subluminal for scalars. 



4. Vacuum Polarization and the Refractive Index 



The propagation of photons at the quantum level is determined by the terms in the 
effective action quadratic in A^(x). This is the vacuum polarization: 

j ^fg^d^x ^/^)d 4 x' A^(x)U^(x,x')A u (x r ) . (4.1) 
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where g(x) = det[—g IM/ (x)]. At the one-loop level, the on-shell condition for the 
photon wavefunction is therefore 

V„F% = -aJ VsMA'n^^x.j'jA^') • (4.2) 

To find the refractive index at one-loop order, we substitute the tree-level form for 
the photon wavefunction inside the integral and take the first term off-shell to give 
an equation for the unknown set of functions oj) in eq. ( J2 . 1 1|) . 

Notice, however, that even the effective action computed to all orders in the 
derivative expansion [28] does not entirely capture the essential physics of high- 
frequency propagation since, as we have shown in ref. [16,17] (see also [6,29]), the 
high-frequency dependence of the refractive index is non-perturbative in the param- 
eter uJ 2 R/m 4i . The analysis of vacuum polarization given here automatically includes 
this crucial non-perturbative behaviour. 



4.1 Vacuum polarization and the Penrose Limit 



The complete one-loop vacuum polarization n^ oop (x,x') receives contributions from 
two Feynman diagrams, as illustrated in Fig. |l|. This gives 

= e 2 g^ u 5^\x — x')Gp(x, x) 



-rrl-loop/ 



+ 2e 2 



D^Gf^x, x')d' l/ Gf(x, x 1 ) — Gf(x, x^d^d^G^x, x) 



(4.3) 



where Gp(x,x') is the Feynman propagator of the massive (scalar) electron. 





Figure 1: The two Feynman diagrams that contribute to the vacuum polarization to order a. 



The Feynman propagator in a general background spacetime can be written in 
the heat-kernel or "proper-time" formalism as 

n i '\ V detA M!/ (a;, x') f°° dT . _ im 2 T+ J_ a i x , x ') n , lA A \ 

GF(X ' X)= Ur^Y" l n(X ' X{Th <4 ' 4) 



subject to the usual m 2 — > m 2 — ie prescription. Here, a(x, x') is the geodesic interval 
between the points x and x': 



a(x,x') = \\ drg^x^x" , (4.5) 
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where x^ = x M (r) is the geodesic joining x = x(0) and x' = x(l). The factor 
det A Aljy (x, x') is the famous Van Vleck-Morette (VVM) determinant, where the ma- 
trix is 

A ( x x >) = g^fe: x ') (46) 

The geometric nature of the VVM matrix and its relation to geodesic deviation is 
explored in detail in Section 5. 

This expression for the propagator has a nice interpretation in the worldline 
formalism, in which the propagator between two points x and x' is determined by a 
sum over worldlines x M (r) that connect x = x(0) and x' = x(T) weighted by expiS[x] 
where the action is 

1 r T 

S[x] = -m 2 T + - dr g^(x)x»x v . (4.7) 
4 Jo 

Here, T is the worldline length of the loop which is an auxiliary parameter that 

X 




Figure 2: The Feynman propagator Gf{x,x') is expressed as a functional integral over paths 
joining x 1 to x. In the limit of weak curvature, R <C m 2 , the functional integral is dominated by a 
stationary phase solution which is the geodesic joining x' and x. 



must be integrated over. The expression ( |4.4|) corresponds to the expansion of the 
resulting functional integral around the stationary phase solution, which is simply the 
classical geodesic that joins x and x' as illustrated in Fig. 0. In particular, the classical 
geodesic has an action S[x] = a(x,x')/(2T) — m 2 T giving the exponential terms in 
( |4.4| ). The VVM determinant comes from integrating over the fluctuations around 
the geodesic to Gaussian order while the term Q(x, x'\T) = 1 + J2^=i a n(x, x')T n 
encodes all the higher non-linear corrections. Notice that these terms are effectively 
an expansion in R/m 2 , so the form for the propagator is useful in the limit of weak 
curvature compared with the Compton wavelength of the electron. Of course, this 
is precisely the limit we are working in here. 



The weak curvature limit R -C m 2 leads to a considerable simplification as we 



now explain. The terms in the exponent in the second term in ( |4.3|) are of the form 

(4.8) 



exp 



im 2 T -( — + — ) v ' - iwV 



For later use, we find it convenient to change variables from T\ and T 2 to T = T% + T 2 
and £ = Ti/T, so < £ < 1. Expressed the other way 

T X = T£, T 2 = T(l-0- (4.9) 

The Jacobian is 

dT x dT 2 r dT f 1 d£ 



Tl If J T*J [e(i-oi 2 ' 

In the limit R m 2 the integral over x' is dominated by a stationary phase 
determined by extremizing the exponent (|4.8Q with respect to x': 



2T^(1 - 



' «9>(x, a;') + w^V" = . (4.11) 



Since 9' M cr(a;, x') is the tangent vector at x' of the geodesic passing through x' and x, 
the stationary phase solution corresponds to a geodesic with tangent vector a <9' M V. 
This means that x and x' must lie on one of the geodesies of the null congruence. If we 
choose x to be the point (u, 0, 0, 0) then x' must have Rosen coordinates (u', 0, 0, 0). 
We call this distinguished null geodesic 7. For these points, it follows that for any 
metric for which dy is a Killing vector, 

dvra(x, x') = u - v! ; (4.12) 

so the V component of ( |4.11| ) becomes 

u' — u , 

+ u = (4.13) 



2Tf (1 - 
and hence 



u' = u - 2u)Ti{\ - . (4.14) 



In the equivalent worldline picture, the stationary phase solution which domi- 
nates in the limit R m 2 describes a situation where the incoming photon decays 
to an electron positron pair at the point v! = u — 2u;T£(l — £) which propagate along 
the null geodesic 7 to the point u and then combine into the photon again, as shown 
in Fig. |3|. This was a key step in the derivation of the refractive index in the worldline 
formalism which we presented in ref. [16,17]. 
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Figure 3: The classical stationary phase solution where the photon travelling along the geodesic 
7 decays to an electron-positron pair at x(u') which both follow the geodesic 7 and then re-combine 
back into the photon at x{u) where u = u' + 2wT'f(l — £). 



In either formalism, the fluctuations around the stationary phase solution are 
governed by the ratio RT which is effectively R/m 2 . In order to set up the expansion 
systematically it is useful to make the following re-scaling of the coordinates 

{u, V, Y a ) — ► (u, TV, VfY a ) . (4.15) 

which implements an overall Weyl scaling while preserving the stationary phase so- 
lution ( |4.14D . After these re-scalings, the geodesic interval Q4.5|) with metric (|2.6|) 
becomes 

a(x,x') = -[ dr \ - 2iiV + C ab (u,0,0)Y a Y b ] + <D{T 2 ) (4.16) 
2 Jo L 

The leading order piece is precisely the Penrose limit around the null geodesic 7 (V = 
Y a = 0). The Penrose limit is the limit of the full metric in a tubular neighbourhood 
of a null geodesic, as illustrated in Fig. |J defined in such a way that it captures the 
tidal forces on the null geodesies that are infinitesimal deformations of 7. (This point 
of view will be described more fully in Section 5.) It follows that to leading order in 




Figure 4: The Penrose limit associated to the null geodesic 7 is the limit of the full metric that 
captures the tidal forces on nearby null geodesies. 

the expansion in R/m 2 , we can replace the metric by its Penrose limit around the 
null geodesic V = Y a = 0: 

ds 2 = -2du dV + C ab {u)dY a dY b , (4.17) 

where C a b(u) = C a b(u, 0, 0). This defines a plane wave in Rosen coordinates. 
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4.2 Geometry of the plane-wave metric 



The fact that the leading-order contribution to the vacuum polarization for an arbi- 
trary curved spacetime depends only on its Penrose limit is a remarkable simplifica- 
tion. As we show below, it allows the derivation of a strikingly elegant expression for 
the full frequency dependence of the refractive index, given purely in terms of the 
VVM matrix. First, we collect some geometrical properties of the plane wave metric 
in both Rosen and Brinkmann coordinates. 



The connection between the Rosen coordinates (u, V, Y a ) and Brinkmann coor- 
dinates (u,v,y t ) involves a zweibein E l a (u), which ensures that the transverse space 
is flat in Brinkmann coordinates. That is, 3 

C ab (u) = E\(u)5 tJ E\(u) , (4.18) 

Then, solving the null geodesic equation in the plane wave metric [16, 17] motivates 
the following coordinate transformation: 



y l = E\Y a , 

dE (4-19) 
v = V + \^E\Y a Y h . 
1 du 



The inverse transformations are therefore 

Y a = y l Ei a 



v = v- ln ijy y 



(4.20) 



where fiy = d u Ei a Ej a plays an important role in the Brinkmann analysis. In partic- 
ular, the zweibein must be chosen in such a way that fiy = Qji. In these coordinates, 
the Penrose limit ( |4.17| ) takes the familiar plane-wave form 

ds 2 = -2du dv - hij {u)y i y j du 2 + dy { dy l , (4-21) 

where the quadratic form is 

hij(u) = R lU3U = - d2E ^ U) E 3 a (u) , (4.22) 

and we have the useful identity h = —d u fl — O 2 . Here, and in the following boldface 
symbols are used to denote 2x2 matrices with Brinkmann transverse i,j indices. 



The Brinkmann coordinates are more fundamental to the distinguished geodesic 
7, Y a = V = 0, than the Rosen coordinates, since they are the geodesic analogues of 

3 Note that the i index on E l a is raised and lowered with 5{j while the a index is raised and 
lowered with C a b(u) and its inverse. Also note that in Rosen coordinates, W g(u) — det_E. 
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Riemann normal coordinates known also as Fermi normal coordinates [24]. In addi- 
tion, the Rosen coordinates are not unique since there are always many inequivalent 
congruences of which 7 is a member. In the following, we find the Rosen coordi- 
nates to be the most efficient for performing the calculation while the final result is 
naturally expressed in terms of the more fundamental Brinkmann coordinates. 

The geodesic interval is particularly simple in Rosen coordinates: 

a(x,x') = (it - u'){V - V) + \A ah {u, u')(Y - Y') a (Y - Y'f (4.23) 

where 



A ab (u,u') = (u-v! 



C-\u")du" 



(4.24) 



A ab are therefore the transverse Rosen components of the full VVM matrix. The 
VVM determinant itself reduces to a determinant over the two-dimensional transverse 
space 

det A^(x, x') = -det A ab {u, u) . (4.25) 



To prove this, we first take r = (u" — u') /(u — u') in the definition ( |4.5|) so that 



a(x,x') = ^^l du"( -2V + C ab Y a Y b ) . (1.26) 



The geodesic equation for the Y a {u) is simply 

d 



, C ab (u)Y b (u) = , (4.27) 
du 



with solution 

- 1 / \ 1 a b 



Y a {u)= [C'\u)] a \ b (4.28) 
for constant ^ b . Integrating this, and using the definition ( [4.24|) , we have 



(Y - Y'f = f U du" [C-\u")] a \ b = (u u>) [A~ 1 (u,u')] ab £ ib . (4.29) 

Ju' 

Hence, the geodesic interval is 



a(x,x') = (u u')(V -V') + ^f{Y- Y'n a 

= _( u _ u ')(y _ V) + lA ab (u,u')(Y - Y') a (Y - Y') 1 



(4.30) 



as claimed. 



In addition, Q(x, x'\T) = 1 in a plane- wave background, which is a manifestation 
of the fact that the propagator is WKB exact. This is entirely consistent with the 
fact that for the original metric the non-leading terms in Q(x, x'\T) are suppressed in 
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the limit m 2 ^> R. The implication of this is that in a general background spacetime 
our analysis is valid in the limits u 3> V~R and m 2 3> R. However, for a plane wave 
spacetime, the results will actually be exact for any R, m and to. 



The eikonal approximation ( 2.10|) for the electromagnetic field A^{x) is similarly 
exact for a plane wave spacetime. Moreover, all the quantities involved have a very 
simple geometric interpretation [16,17]. Specifically, the amplitude is 



A(x) = (detE ia (u)) 1/2 
and the non-vanishing components of the polarization vector e^w are 



i(i)a{x) = E ia {u) 



The tree level contribution to the mass shell condition 
for small is then simply 



at a point x 



(4.31) 

(4.32) 
{u, 0,0,0), 

(4.33) 



4.3 Refractive Index 



We now complete the calculation of the vacuum polarization and refractive index, 
working from here onwards in a plane wave background. Returning to the expression 
( |1.3|) for the vacuum polarization, we find that the first Feynman diagram in Fig. [l] 
gives the following contribution to (|4.2f ): 



-A(u)E ia (u) / — le 
7T Jo T 2 



■im 2 T 



(4.34) 



By itself, this contribution is divergent but we shall find that it cancels a divergence 
in the second term. 



The contribution to the on-shell condition ( |4.2| ) from the second Feynman dia- 
gram in Fig. [I] is then 



8e 2 / ^/g(x')d 4 x' d„G F (x, x')& v G F (x, x') — Gp{x, x')d^d' u GF(x, x') 



v -iuV(x') 

fc (i) e 

(4.35) 

In Rosen coordinates, we take x = (u, 0,0,0) and x' = (u' ,V ,Y' a ). What remains 
is to integrate over x'^ . The integral over V is trivial and leads to a delta function 
constraint 



dV' exp 



(u' - u)V 
L2iT£(l-0 



4ttT^(1 - £)6(u' -u + 2uT£(l - 0) (4.36) 
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which saturates the v! integral. This simply enforces the condition ( 4.14 ), the sta- 
tionary phase solution becoming exact for the plane wave background. 



Since only has non-vanishing components in the Y a directions and the in- 
tegrals over the Y' a are Gaussian, it follows that ( 4.35|) is only non-vanishing if the 
derivatives <9 M lie in the directions d a . Using this fact, the Y' a integrals are of the 
form 

d 2y> 9 ^Y'-A(u,u')-Y> _d_ 4T( \_ n Y'-A(u,u')Y' _ 7r£(l - Q A ab (u, u') 

(4.37) 

This is where the advantage of performing the calculation in Rosen coordinates is 
clearest, since these coordinates automatically exhibit the simple form (|4.23|) for 
the transverse sector of the geodesic interval. (The corresponding expression in 
Brinkmann coordinates is given in Section 5.) Noting that 



A ab (u, u')Ei{u') = E 3 a (u)A ij (u ) u') 



(4.38) 



and the fact that 



'd et A ab (u,u') 
9{u) 



{&etE ia (u')) 1/2 = (detE ia (u)) v Vdet E^(u)A ab (u, u^E^u' 



1/2 



A(u) J det Ay (m, u' 



(4.39) 



we find the contribution to the mass shell condition ( [4.2[ ) is 

dT rl 



-A{u)E\{u 



7 , 2 , diie~ m T A ij (u,u') x /detA ij (u,u' 



u'=u-2u)T £,{!-£,) 



(4.40) 



Summing over the two contributions to ( [4.3| ) gives the complete one-loop term 
in ( [4.2|) and since the tree- level contribution is ( |4.33| ), we can extract the matrix of 
refractive indices: 



a 



n[u\ oj 



27VL0 2 



dT 

~J>2 



ie 



-im z T 



X / di 
'0 



1 - A(u, u) ^/detA(u,M') 



(4.41) 



J u'=u-2u>T£(l-£) 



Here, and in the following, A represents the 2x2 matrix Ay with Brinkmann 
coordinate indices. 



This is our principal result for the refractive index. It is remarkable that the full 
frequency dependence of the refractive index/phase velocity for photons propagating 
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in an arbitrary background spacetime can be expressed in such a simple and elegant 
way. The key insight, that the quantum effects on photon propagation are determined 
by the geometry of geodesic fluctuations around the classical null trajectory and are 
therefore entirely encoded in the plane-wave Penrose limit of the original spacetime, 
explains why the final result should depend so simply on the VVM matrix only. 

Notice that the result for the refractive index at a point x = x(u) only depends 
upon data associated to the classical null geodesic x(u — t), < t < oo, i.e. on the 
portion in the past relative to x. We can write 



a 



2nuj 







(4.42) 



with 



F(u: z) 



dt 
1? 



le 



—izt 



1 - A (w, u - t) J det A (u, u-t) 



(4.43) 



where we changed variables from T to t = 2a;£(i — £)T. The result we have obtained 
is strictly valid for uo real and positive. Also notice that the definition ( ^4.43| ) has the 
form of a Fourier Transform of a function which vanishes for t < 0. 



X 



X 



X 

X 
X 

X 



XXX 

■ — y~ 



Figure 5: The integration contour in the i-plane that defines the physical values of the refractive 
index for real positive lj. Crosses represent branch-point or pole singularities which generically lie 
on the real axis but in some examples lie also on the imaginary axis. 



We will show in Section 6 that as a consequence of singularities in A corre- 
sponding to conjugate points on the null congruence, the integrand has branch-point 
singularities on the real t axis and so the t integral must be defined by some prescrip- 
tion. The ie prescription that we have chosen in ( [4.43Q ensures that the refractive 
index becomes trivial in the flat-space limit. The t integration contour is illustrated 
in Fig. m 
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The low frequency behaviour of the refractive index follows readily by expanding 
the VVM matrix in powers of t, since the effective expansion parameter is u 2 R/m 4 . 
Expanding 

oo 

A(u,u-t) -> l + ^s (n) (w)t 2n , (4.44) 

n=l 

where s^ n \u) oc R n , we identify the term linear in the curvature as 

4V) = -\Rmuj(u) . (4.45) 
Substituting this expansion into ( [4.42 ), we therefore find to leading order 

mj (u; u) = 5ij + 36Q " m2 (RuJij - 2R uluj ) + ^o(^J^j . (4.46) 



in agreement with the result (3.10) derived from the effective action. 



5. Geodesic Deviation and the VVM Determinant 



The Van Vleck-Morette determinant plays a central role in determining the refractive 
index and its analytic structure. This is because the VVM matrix controls the 
geometry of geodesic deviation. Since this is such an important part of our analysis, 
in this Section we present a detailed account of this geometry, mostly from the 
viewpoint of Brinkmann coordinates. 

We start with the definition. Fix two points x' and x in spacetime and consider 
the following functional integral 

j \dx v^MK^ 1 , (5.1) 

where the action is 

S[x] = \ [ drg, v {x)x»x v (5.2) 
Jo 

and the function x{r) has boundary conditions x(0) = x' and x(l) = x. The WM 
determinant arises from integrating the fluctuations to Gaussian order around a 
stationary phase solution of the equations of motion. These are precisely the geodesic 
equations. If we denote by x(r) a (usually unique) geodesic that passes through 
x 1 = x(0) and x = x(l), the equations for the fluctuations are 



W v 6x v (t) = , 



(5.3) 



where V^ v is the second order differential operator 

V ^ = <^TT2 + R^uxx a i x . (5.4) 

Here, D / Dt are absolute derivatives alongx(r), i.e. Dx^{r) / Dr = d T x^+T^{x) va x v x a '. 
The VVM determinant is defined as the functional determinant detX"V It is eval- 
uated directly in [35] (see also [36]) by discretizing the functional integral and then 
taking a continuum limit to yield the finite determinant 

detA^(l) , (5.5) 

where A fJ, u (r) is a solution of the Jacobi equation 

V» v A\{t) = , (5.6) 

subject to the boundary conditions 

A"„(0) = , = 8\ . (5.7) 

To understand this in more detail and connect to the previous definition of the 
VVM determinant, we now specialize to Brinkmann coordinates and consider the 
tidal forces on null geodesies that are infinitesimal deformations of the distinguished 
null geodesic 7. 4 These nearby geodesies are described by the "Jacobi fields" in the 
neighbourhood of 7 which, given their interpretation as Fermi null coordinates, can 
simply be identified as the transverse Brinkmann coordinates y l (u) of null geodesies 
in the plane wave metric ( |4.21|) . Their evolution is described by the geodesic deviation 



equation, specializing ( |5.3| ) and choosing u as the afline parameter: 

d 2 y l (u) 



dv? 



-R\ JU (u)y\u) . (5.8) 



The solution of eq. (|5.8|) determines the coordinates y(u) in terms of initial data 



y{u') and y(u') at a fixed point u', i.e. 

y\u) = B l 3 {uy)yi{u') + A^(u,u')y^u f ) . (5.9) 

It follows immediately that the matrix functions A(u, u') and B(u, u'), with elements 
Aij(u, u') and Bij(u,u'), satisfy the geodesic deviation equations 

A + hA = 0, B + hB = 0, (5.10) 



4 It is important to realize that these nearby geodesies do not necessarily lift to geodesies of the 
full metric: this is a global issue of integrability that is irrelevant to our discussion. 
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where h has elements hij(u) = Ri U j U , with boundary conditions A(u',u') = 0, 
d u A(u, u')\ u=u i = 1, B{u' ,u') = 1 and d u B(u, u')\ u=u i = 0. Using the zweibein, 
it is easy to see that these functions satisfy the consistency relation 

B l3 + A lk n k 3 {u') = E ia {u)E j a (u') . (5.11) 



Two special choices of boundary conditions for y l (u) are of particular interest: 

(i) y l (u) = y i (u) = A l 1 (u,u')y\u') , (5.12) 

with, as always, = and d u Aij = 5ij at u — u' . This describes a "spray" of 
geodesies [26] passing through a point y{u') = and determines the function Aij 
which, as we show below, is related very simply to the inverse of the VVM matrix 
Aj,j. In addition, as we prove below, A(u,u') has the anti-symmetric property 

A(u,u') = -A(?/,u) T . (5.13) 



(ii) y\u') = y\u) = B^(u,u')y^u') , (5.14) 

with Bij = 5y and d u B^ = at u = u' . This is the choice [25] appropriate to a 
geodesic congruence with neighbouring geodesies parallel at u' . 

The geodesic deviation functions A and B determine the geodesic interval for 
the plane wave in Brinkmann coordinates. Analogous to the Rosen expression ( |4.2 
we have 



u' 

— iu — u) I (V 



(5.15) 



tr(x, x') = I ^ du" (2v + h l3 (u")yY - y%^j 

using the geodesic equation y % + h l jy^ = 0. Substituting ( |5.9|) now gives 

a(x, x') = -(u- v!) {{v - v') + y\u)Ajl{u, u')y j (u') 

+ ly\u)A^(u',u)B k ,(u',u)y^u) - \y\u')A^{u, u')B k j(u, u')y j {u') 

(5.16) 

The transverse Brinkmann components Ay of the VVM matrix, defined by 

A ti (u,u r ) = X '} r (5.17) 
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and therefore 

A{u, u') ={u- v!) (A-^u, u')) T . (5.18) 

Yet another interpretation [26] of the VVM determinant is as the Jacobian for 
the change of variables between specifying a geodesic by giving two points — x(u') 
and x(u) — through which it passes and giving one point and the tangent vector at 
that point: x{u') and x(u'). Since we can think of k u (x,x') = da(x,x')/dx w as 
the tangent vector at x' of the geodesic that goes through the two points x' and x, 
normalized so that the affine parameter between x' and x goes from to 1, we see 
from flyD that 

A/(x,x')=d fl k v (x,x') . (5.19) 

With this normalization, 

k\u,u') = {u-u')y\u') (5.20) 
and so from ( |5.12|) we have the Jacobian matrix 

d -M^l = ( u-u-)AJ,Hu,u'). (5.21) 



The equivalence of the Rosen and Brinkmann expressions ( |4.23| ) and ( |5.16| ) for 
the geodesic interval is readily established once the following identity is proved: 



Aij{v.,v.') = E ia {u') I du"[C-\u")] ab E jb (u) . (5.22) 

Ju> 

The proof is as follows. Notice that the zweibein Ei a (u) is a particular solution of 
the geodesic equation (fj.iop . A{u,u') also solves this equation and so it follows that 



E T A -E A = K , (5.23) 

where K is a constant matrix. Using the fact that f2 = EE~ X is a symmetric matrix 
allows us to write 

^(.Er 1 ^) = E-^E-^K = C l K , (5.24) 

where C(u) is the non-trivial part of the metric in Rosen coordinates ( |4 . 1 7| ) . In- 
tegrating this equation and imposing the boundary conditions A(u', u') = and 
d u A(u,u')\ u=u > = 1, gives 

A(u,u') = E(u) [ du"C-\u")E(u'Y . (5.25) 

Ju' 
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which in components is ( 5.22 ). Notice that the symmetry ( 5.13 ) is manifest. A 
similar construction with the alternative boundary conditions determines B(u,u') in 
the form ( |5.11| ). 5 



Finally we relate these results to the optical scalars in the Raychoudhuri equa- 
tions which describe the geodesic flow. It is convenient to start from an alternative, 
but entirely equivalent, description of geodesic deviation. In this approach, the evo- 
lution of the Jacobi fields y l (u) is determined by requiring that their Lie derivative 
CkV 1 vanishes along the geodesic with tangent vector i.e. 

C k y l = k-Vy l ~ (V^V = . (5.26) 

This implies the parallel transport equation 

d u f = tfrf (5.27) 

where Qij = Vjki. Notice that since the geodesic tangent vector in Brinkmann 
coordinates is /c M = (1, (|f2 + QL 2 )ijy l yi , fi^-y 3 '), this is consistent with the original 
definition Qij = Ei a Ej a . It then follows from ( |5.9|) that 

du^B^u^ + A^u')^^')) =fi ii («), (5.28) 

which is clearly consistent with ( p. lip . 



The matrix is the fundamental object from which the optical scalars are 
defined. We have 6 

flij(u) = \9{u)8ij + &ij(u) + u>ij(u) , (5.29) 

defining the expansion 8, the shear <Jij and the twist il>ij. The corresponding scalars 
are a 2 = ft^ft^ — \6 2 and Cj 2 = Vt^^Vt^ . The twist vanishes in all cases considered 
here, so Qij is symmetric. Eq. (|5.28 ) therefore implies: 

B ij (u,u')+A ik (u,u , )n k j (u') = exp / du" (±96$ + . (5.30) 

Ju> 



Note that the optical scalars depend on the choice of boundary conditions im- 
posed on y l {u). A particularly relevant choice is the "geodesic spray" condition 
considered above. In this case, ( |5.28| ) simplifies to 



cUogA = Q(u) 



(5.31) 



3 Also notice that if we were to evaluate the vacuum polarization directly using the Brinkmann 



expression for a(x, x 1 ) rather than the simpler Rosen form, ( 5.11 ) is essential in simplifying the 
transverse integrals and ensuring that the elegant Rosen result ( 4.37 ) is reproduced. 

6 Here, we follow the conventions of Wald [27]. There are therefore some factors of 2 different 
from the Chandrasekhar [37] conventions used in refs. [16,17]. 



- 22 - 



Taking the trace gives the important identity 

d u log det A(u, u') = 9(u, u') , (5.32) 

where we display the u' dependence on the r.h.s. explicitly as a reminder of the choice 
of boundary condition, just as in the notation k l {u, u') in (|5.21| ) for the tangent vector. 




Figure 6: Conjugate points x{u) and x(v!) on the geodesic 7 through which infinitesimal defor- 
mations of 7 (or finite deformations in the Penrose limit) also pass. 

In general, if there are two points x(u) and x(u') on a geodesic 7 for which there 
exists a family of geodesies infinitesimally close to 7 which also pass through x(u) 
and x(u'), then these are said to be conjugate points. As we now show, conjugate 
points play a crucial role in determining the analyticity properties of the refractive 
index. For the plane wave, conjugate points correspond to solutions of the geodesic 
equation with y l (u) = y l (u') = 0. It follows from the discussion above that this 
implies det A(u, u') = 0. In turn, this implies that at these points, the VVM de- 
terminant detA(-u,-u') has a singularity. This establishes a direct link between the 
analyticity structure of the refractive index ( |4.41| ) and the geometry of conjugate 
points. Moreover, this geometry is entirely encoded in the geodesic deviation matrix 
A(u,u') or equivalently the VVM matrix A(u,u'). 




Figure 7: The expansion 9(u, u 1 ) describes the rate-of-change of the proper area defined by the 
spray of null geodesies that pass through x(u'). The null energy condition implies that the expansion 
monotonically decreases, except at conjugate points. 
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A particularly important observation is that the existence of conjugate points is 
generic (see e.g. [27]). This follows from the Raychoudhuri equations: 



dj = -\9 2 - a^a 11 - R uu ^ 
As a consequence of the null energy condition, R uu > 0, ( |5.33|) implies the inequality 



dj(u, u') + |0(u, u'f < . (5.34) 

The significance of this is that d u 9(u,u') < so that 9(u,u') generally decreases 
monotonically with u. (Of course, this is violated at the singularities where 9(u,u') 
jumps from — oo to oo.) If at some point u = u, 9(u,u') is negative, say — |A|, then 
inevitably 9{u,u') — > — oo at some finite u < tt + 2/|A|. The proof is simple. In order 
to attempt to avoid the singularity \d u 9(u, u r )\ should be as small as possible. In 
other words, we should saturate the inequality ( |5.34| ), with the solution 



9(u,u') = . (5.35) 

u — u — 2/ \ A\ 

Hence, there must be a conjugate point at some x(u) with u < u + 2/|A|. At the 
conjugate point, 9{u,u') jumps discontinuously from — oo to oo and then begins its 
descent again. Notice that as u — ► oo the expansion must go asymptotically to zero. 




Finally, notice that by diagonalizing the shear tensor and combining the two 
Raychoudhuri equations (|5.33| ), we can characterize the null congruences by whether 
the geodesies focus in both transverse directions (specified by the eigenvectors of <7jj) 
or have one direction focusing and one defocusing. (These were labelled as Type I 
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and Type II respectively in refs. [16,17].) As shown above, the null energy condition 
prohibits the existence of a third case with defocusing/defocusing. The focusing 
directions give rise to conjugate points and corresponding singularities of the VVM 
determinant A.(u,u') on the real axis; defocusing directions, on the other hand, are 
associated with singularities on the imaginary axis, as illustrated in the example of 
general symmetric plane waves in Section 7.2. 



6. Analyticity and Causality 



As noted at the end of Section 4, the singularities in the VVM determinant induced 
by the existence of conjugate points gives rise to a novel analytic structure for the 
refractive index in the complex u plane. This is a generic effect which will also 
affect more general scattering amplitudes. As we shall see, it means that in curved 
spacetime some of the conventional axioms and assumptions of S'-matrix theory and 
dispersion relations need to be re-evaluated, with far-reaching physical implications. 



6.1 Analytic structure of the refractive index 



Returning now to the expression ( 4.42 ) for the refractive index in terms of the VVM 
matrix, 



a 



n[u] lu 



J dee(i-o^(; 



m 



it is clear that the t- integral defining J~{u\ z) in ( 4.43) has branch-point singularities 
on the real t-axis whenever x(u) and x(u — t) are conjugate points. After integration 
over t, these singularities give rise to cuts in 7-{u\ z) from to oo in the complex 
z-plane. The refractive index n(u; lu) is therefore a multi-valued function of uj with 
branch points at and oo. 



The presence of branch-point singularities on the real t-axis means we have to 
give a prescription for the contour of the t- integration in ( |4.43| ). For z real and 
positive, as in Section 4, we define: 



zeR>o 



dt 



-izt 



1 - A (u, u - t) ^detA(u,u - t) 



(6.2) 



with the contour as illustrated in Fig. |5|. With this choice, the t-integral can be per- 
formed by rotating to the negative imaginary t-axis where the integral is convergent 
due to the damping of e~ tzt as t — > —too for Re z > 0. Moreover, since the integral 
is over Re t > 0, it receives support only from that part of the null geodesic to the 
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past of x(u), i.e. from A.(u,u — t) with Re t > 0. Intuitively this is what one would 
expect for in a causal theory. 



Similarly, for z real and negative we should define 



T-{u\ z) 



2eR<o 



f-°°+ ie dt r / 

/ — ie~ lzt 1- A(u,u-t)JdetA(u,u-t) 

Jo t L 



(6.3) 



This contour avoids the singularities on the negative t-axis, while rotation of the 
contour towards the positive imaginary t-axis leads to a convergent integral for Re 
z < 0. This time, the t-integral has support only from the section of the null geodesic 
in the future of x(u). Again, this is as required for causality with Re a; < and is 
consistent with the usual flat-space limit. 

The next step is to specify the physical sheet for the multi-valued function 
ZF(u; z) defining the refractive index. First, we choose to run the branch cuts in 
the complex z-plane from to ±oo just above (below) the positive (negative) real 
axis respectively. 7 The physical J-(u\ z) is defined as the analytic continuation of 
J- '+(«; z) from real, positive z into the lower-half plane and of ZF-(u; z) from real, 
negative z into the upper-half plane. That is, 

, - _ , — 7r < arg z < , . 

T{u-z) = { J) ; - (6.4) 

< arg z < 7r . 




This implies the corresponding analytic structure for the refractive index n(u; u) 
itself in the complex tu-plane, illustrated in Fig. [{J Since z is essentially the inverse of 
u, the upper-half plane in z maps into the lower-half plane in uj and vice-versa. The 
physical refractive index is therefore given by the analytic continuation of n + (u; un- 
defined using !F + {u\ z) — into the upper-half plane and of n_(u; uj) into the lower-half 
plane 

There may also be further singularities in ^F + (u;z) or F-(u;z) individually 
(e.g. we will find examples in the next section where J~+(u; z) has poles on z G R < 0) 
but these lie off the physical sheet defining T(u\ z) itself. 



Across the cuts, there will be a discontinuity in J~(u; z) or n(u; uj) and we define 
Disc J-(u\ z) and Disc n(u; uj) as the discontinuities across the appropriate cuts taken 
in the anti-clockwise sense. These discontinuities play an important role in dispersion 
relations. 

7 This applies also to the cuts arising from any branch-point singularities occurring off the real 
t-axis, for example the singularities on the imaginary i-axis in the general symmetric plane wave 
example discussed in Section 7. 
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CO 



Figure 9: The physical sheet for n(u;uj) defined by analytic continuation from the real positive 
axis into the upper-half plane and from the real negative axis into the lower-half plane. Since z 
is inversely related to u> the upper-half of the u> plane maps to the lower-half of the z plane, and 
vice- versa. 



In the simplest case, where there are no singularities in the complex t-plane apart 
from those on the real axis (as realized in the conformally flat symmetric plane wave 
example discussed in Section 7.1), we can evaluate Disc J-{u\ z) across the cut along 
z G R > by rotating the contour defining 7-^{u]z) to wrap around the positive 
t-axis as shown in Fig. |H]. That is, 





Figure 10: Wrapping the contours C_ and C+ defining T '-(u; w) and F+(u; lo) around the branch 
point singularities on the positive real t-axis. The integral with the resulting contour C gives the 
discontinuity of J~(u; tv) across the cut along the positive real w-axis. 



Disc J-{u\ z) 



zeR>o 



T-{u\ z) — J 7 + (u; z) 



dt 



-ie 



dt 
72 



-izt 



dt 
72 



ie~ lzt [..] 



(6.5) 



1 - A(«, u - t) J det A [u, u-t) 
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Indeed, for the conformally symmetric plane wave background, the singularities on 
the real i-axis are actually poles, so Disc J-{u\ z) can be evaluated from the contour 
C simply as the sum of the residues. This example is worked out explicitly in Section 
7.1. 

An important special case arises when the background is translation invariant 
with respect to the coordinate u along the null geodesic. Since the WM matrix is 
symmetric in its two arguments, we then have 

A(u,u-t) = A(u-t,u) = A(u,u + t) , (6.6) 

so the integrand [. . .] in tF{z) (which is of course then independent of u) is an even 
function of t. In turn, this implies 

J- -{z) = —J- + (—z) (translation invariance) (6.7) 

which is a result of special significance in the analysis of dispersion relations. 



6.2 Kramers-Kronig dispersion relation 



The Kramers-Kronig dispersion relation is an identity satisfied by the refractive 
index, or vacuum polarization, in QED. 8 Its derivation depends critically on the 
analyticity properties of the refractive index and shows in a simple context the sort 
of changes to conventional ^-matrix relations and dispersion relations which will 
occur due to the novel analytic structure of amplitudes in curved spacetime. 



To derive the Kramers-Kronig relation, we integrate n(u;uo)/uo around the con- 
tour shown in Fig. 11, As explained in the following section, causality imposes two 
fundamental properties of the refractive index: (i) n(u; uo) is analytic in the upper- 
half cu-plane, and (ii) n(u; uo) is bounded at infinity. 9 Assuming these properties, we 
have 

duo . . f duo _ f°° duo . . 

— n{u]uo) + / — n(u]uo) + V \ — n{u;u) = 0. (6.8) 

U JCo U J-OD ^ 

which implies 

n(u; 0) — n(u; oo) = ^- V [ —n(u;uo). (6.9) 



m .7_oo uo 



8 Contrary to some recent claims in the literature [11], the Kramers-Kronig relation is equally 
valid in relativistic quantum field theory as it is in non-relativistic settings; for example the proof 
in QFT is presented in Weinberg's textbook [15]. 

9 Notice that this is weaker than the condition that n{u;uS) — * 1 as \uj\ — * oo. 
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Figure 11: The contour in the complex w-plane used in the derivation of the Kramers-Kronig 
relation for n(u\ uj). The second figure shows the equivalent contour in the complex z-plane relevant 
for ^(u; z). 

Provided the causality properties are satisfied, the Kramers-Kronig relation in the 
form ( p.9|) is always valid. Note that the principal part integral is over a contour 
lying just above the real axis, so can be written as 




(6.10) 



Now, in the conventional flat-spacetime derivation [15], translation invariance 
implies that n{uj) is an even function of uj. This implies that n(— uj + ie) = n(uj — ie), 
so the r.h.s. of ( |6.1Up becomes the discontinuity of n(uj) on the positive real cu-axis: 



n(0) — n(oo) = — [ — Discn(co>) (translation invariance) . (6.11) 

m J uj 

Finally, in flat-spacetime QED, the refractive index satisfies the property of real 
analyticity, n(uj*) = n(uj)*. This is a special case of the basic S'-matrix property of 
hermitian analyticity [38] which is satisfied by more general scattering amplitudes. 
With this assumption, we can replace the discontinuity in (|6.11| ) by the imaginary 
part of the refractive index, since for real uj it implies n(uj — ie) = n(uj + ie)*, leaving 

2 f°° duj 

n(0) — n(oo) = — — Imn(u) (real analyticity) . (6.12) 

Jo v 

This is the standard form of the Kramers-Kronig relation. Since the optical theorem 
relates the imaginary part of forward scattering amplitudes to the total cross section, 
the r.h.s. of ( 6.12j ) is positive under conventional QFT conditions. This would imply 
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n(0) > n(oo), consistent with a subluminal low-frequency phase velocity, which is 
the usual dispersive situation. 10 



In curved spacetime, however, the assumptions leading to the second (|6.11|) and 



third ( |6.12| ) forms of the Kramers- Kronig relation need to be reassessed. We still 
maintain the causality conditions, that n(u; uj) is analytic in the upper-half tu-plane 
and bounded at infinity, so the primitive identity (|6.9|) is always satisfied. 



Now, in our case, because of the cuts along the real cu-axis and the definition of 
the physical sheet, the r.h.s. of ( |6.10| ) actually involves the function n + (u; uj). Then, 
if we are in a special case where we have translation invariance along the geodesic, 
so that the VVM matrix A.(u,u') is a function only of (u — u'), we have from ( |6.7| ) 
that 

n-{yj) = uj) , (6.13) 

taking into account the extra 1/uj factor in front of the integral (|6.1| ) for the refractive 
index in terms of F{z). The r.h.s. of ( |6.10| ) therefore involves 

n + {uj + ie) — n + (— uj + ie) = n + {uj + ie) — n_{uj — ie) 

TV < \ ( 6 - 14 ^ 

= Disc n(uj) . 

So translation invariance in u would imply that the second form ( |6.11| ) of the Kramers- 
Kronig relation holds even in curved spacetime. 

However, as we shall see in a number of examples, it appears that real analyticity 
of n(u; uj) is lost for QED in curved spacetime. This stems from the need to define the 
physical sheet for n(u; uj) as in Fig. |9| in terms of both n + (u; uj) and n_(w; uj), with 
the cuts in the complex cu-plane originating directly from the geometry of geodesic 
deviation and the VVM matrix. So the third form (|6.12|) of the Kramers-Kronig 



relation does not hold in curved spacetime. Notice, however, that this does not 
imply there is anything wrong with causality or micro causality. 



6.3 Causality and the refractive index 



In this section, we show that the two conditions on the refractive index assumed 
in the derivation of the Kramers-Kronig relation, viz (i) n(u; uj) is analytic in the 
upper-half cu-plane, and (ii) n(u;uj) — > 1 for large \u\, are necessary conditions for 
causality. The first is a consequence of requiring that the refractive index at x(u) 
only depends on influences in the past light cone of x(u); the second imposes the 

10 For examples in atomic physics where the system exhibits gain and Imn(w) is negative, see 
ref. [8]. See also Sections 6.4 and 7. 
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condition that the wavefront velocity (which has been identified in previous work as 
the relevant speed of light for causality) is c. 



There are many ways to see how the connection between analyticity and causality 
arises, both at the level of the refractive index and more generally in the construction 
of Green functions obeying micro-causality. The essential technical feature is the 
theorem that the Fourier transform f(z) of a function f(t) which vanishes for t > 
is analytic in the upper-half complex z-plane. This argument naturally appears in 
some guise in all the discussions linking causality with analyticity. 




Figure 12: A sharp- fronted wavepacket propagating with respect to the light cone. 

An illuminating illustration is to consider the propagation of a sharp-fronted 
wave packet. Consider the one-loop corrected modes (|2.10| ). Suppose that in the 
distant past, u — > — oo, we build a sharp-fronted wave-packet propagating in the u 
direction, by taking a Fourier Transform of the modes ( p. 1 0| ) : 11 



A» = I dwZ{u)e^{u)e-^ v -^ uvj)) , (6.15) 

with 



/u 
du' [n(u r ; u) - l] (6.16) 
-oo 

In the limit, u — > — oo, the refractive index must approach 1 and so the condition 
that (|6.15D be sharp-fronted, that is vanishing for V < 0, is that Z[lS) is analytic 



in the upper-half plane. This is because when V < we can deform the integration 
contour in ( |6.15| ) from the real axis into the upper-half plane and out to the semi-circle 
at infinity on which the integrand vanishes. 12 This contour deformation argument 
manifests the usual link between causality and analyticity. It follows that at finite 



11 Implicitly we have been assuming that in the distant past the curvature is turned off and the 
refractive index is asymptotically 1. Later we will be able to remove this restriction when we discuss 
Green functions. 

12 In order for the wavepacket to be properly defined in the presence of the one-loop correction, 
the integration contour must lie just above the real axis to avoid any non-analyticities of n(u;u>) 
on the real axis. 
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u, the wavepacket will remain sharp-fronted and vanishing outside the light cone 
V < provided that (i) the refractive index n(u;u), is an analytic function of uj in 
the upper-half plane; and (ii) n(u; uj) — > 1 as \u\ — > oo. If both these conditions are 
satisfied then the sharp- fronted disturbance will propagate causally. 13 



To see how these conditions are realized here, consider the expression ([4.42), (4.43) 



for the refractive index in terms of an integral over t of a function of the VVM matrix 
A.(u,u — t). Causality requires this to depend only on the part of the geodesic to the 
past of x(u), which is guaranteed by the integral being only over t > 0. But since 
J- + {u\ z) then has the form of a Fourier transform of a function which vanishes for 
t < 0, it follows from the above theorem that ^F+(u; z) is analytic in the lower-half 
z-plane. In turn, this implies analyticity of n(u; uj) in the upper-half cu-plane. 

For the second condition (ii), note that the wavefront velocity is the high- 
frequency limit of the phase velocity. 14 For large |u;|, we have 

n(u-uj) = l —T(u;0) + --- (6.17) 

12-kuj 

So a sufficient condition for the second requirement to be satisfied is that the integral 
J-{u\ 0) is finite. 15 In particular, this requires that the integral 



T+{u- 0) 



dt 



1 — A (it, u — t) J det A (u, u-t) 



(6.18) 



is convergent. This is actually guaranteed by the fact that, as we have already 
mentioned, implicitly we are assuming that the space becomes flat in the infinite past 
and future in order that the one-loop corrected photon modes be defined consistently 
in the whole of spacetime. In that case, for large t, A(w, u — t) — ■> 1 as t — ► oo. 16 



13 Notice that causality might be respected if n(u; uj) approached a constant for large \lu\ but with 
respect to a modified light cone given by V = i?(u; oo). This would rely on the space being suitably 
"causally stable" [6, 10,25]. The examples that we find do not have this property and so we will 
not pursue this idea. 

14 This is proved in ref. [12] (see also [6]) for a very general class of wave equation; this proof may 
not, however, be sufficiently general to cover the full vacuum-polarization induced wave equation 

15 It might have been possible for !F{u; z) to have a simple pole at z = 0, in which case the high 
frequency phase velocity is finite but different from c. But as we have already mentioned this does 
not occur. 

16 We can also discuss spaces which do not become flat in the infinite past and future. In that 
case, the relevant problem to consider is an initial value problem and this inevitably involves the 
Green functions, a topic that we turn to in Section 9. 
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6.4 Dispersion and Imn(w) 



In general, the refractive index is a complex quantity. While the real part determines 
the local phase velocity at a point, the imaginary part describes dispersion. To see 
this we note that the probability density of the photon wavefunction, with the one- 
loop correction include, is 



9(u) 



-1/2 



exp 



u / du' lmriij(u';uj) 



(6.19) 



The pre-factor here is just the volume effect one would expect in curved spacetime. 
The exponential term, on the contrary, determines the dispersive effect of space- 
time. In general, we would expect that the eigenvalues of the imaginary part of the 
refractive index would be > so that spacetime would act like an ordinary disper- 
sive medium with the total number of photons being depleted as they propagate by 
conversion into real e + e~ pairs. 



However, when we look at some of the examples that we have considered, in 
particular the Ricci-flat symmetric plane wave whose refractive index is plotted nu- 
merically in Fig. [15| and the weak gravitational wave in Fig. [16|, we see that the 
imaginary part of the refractive index can be negative. Indeed, in the gravitational 
wave example it oscillates sinusoidally with the frequency of the background wave. 
Apparently, in these cases, spacetime acts as an amplifying medium for photon prop- 
agation. 



In many ways, this effect is similar to that studied in some atomic physics exam- 
ples in [8] . It was shown there that for certain three-state A-systems interacting with 
coupling and probe lasers, the refractive index can be arranged to be of the form 

n(u) = 1 , (6.20) 

uj — ujq + ry 

where c^o is a characteristic frequency of the coupled laser-atom system. In the usual 
dispersive case, we would have A > so that Ren(0) > 1 and the low-frequency 
propagation is subluminal while Imn(w) > 0. However, in Raman gain systems we 
can arrange to have A < 0, resulting in a superluminal Ren(0) < 1 with Imn(w) < 0. 
The negative imaginary part indicates that the probe laser is amplified (taking energy 
from the coupling laser), with the system acting as an optical medium exhibiting gain. 

For our purposes, the important point is that in this model, n(u) is characterized 
by a simple pole at uj = ujq — ij. This must be in the negative imaginary half-plane 
to be consistent with causality. It then follows that the existence of Imn(w) < 
is necessarily linked to superluminal low-frequency propagation. In the examples 
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below, we also find that the occurrence of an imaginary part for the refractive index 
is correlated with the occurrence of singularities, in this case branch points, in n(u>) 
off the real axis but in the causally-safe half-plane. In turn, the location of these 
singularities is intimately related to the location of singularities of the VVM matrix in 
the complex w-plane, with polarizations exhibiting Imn(w) < and a superluminal 
Ren(O) < 1 corresponding to the diverging direction of the null geodesic congruence. 
However, we should be cautious about over-interpreting our results in this way, since 
the actual QFT results for the refractive index are significantly more complicated 
than (lOoQ . 

It is also important to recognize that this amplification occurs for photons of 
high frequency ~ 1, i.e. u ~ j-y- where L is the curvature scale. It is not a 
long-range, infra-red effect with photon wavelengths comparable to the curvature, 
ijj ~ -|. Rather, the effect seems to be a kind of emission of photons induced by 
the interaction of the incident wave with quantum loops in the curved spacetime 
background. However, the details of this mechanism remain to be fully understood. 



7. Example 1: symmetric plane waves 



To illustrate these general results, we now consider some simple examples. The 
simplest case is when the background has a Penrose limit which is a symmetric plane 
wave. In this case, the matrix functions R U i U j(u) are independent of u. They can 
immediately be diagonalized, R U i U j = —a 2 5ij, with the <7j constant. The metric in 
Brinkmann coordinates therefore takes the form 

ds 2 = -2du dv - a 2 y { y i du 2 + dy i dy { . (7.1) 

The curvatures are R uu = J^of an d Ruiui = of- 



The VVM matrix can be determined by solving the Jacobi equations (|5.8|) and 



([5.101). This gives y l {u) = C\ cos(o"j-u + c-i) and implementing the appropriate bound- 



ary condition in (|5.10| ) selects 

sin<7j(w — it') 



Hence 



Aijiu, u') = 5ij ^ '- . (7.2) 



Aijiu, u ) = 5 i:j - . (7.3) 

smav [u — u') 
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The matrix of refractive indices is independent of u and diagonal with elements given 
by eqs. - Q with 



dt 



ie 



-izt 



Ojt 

sin cr.jt 



n 

3=1 



sin <jjt 



(7.4) 



and similarly for !F-{z). 17 Notice that in general the integrand has branch-point 
singularities on the real t axis since at least one of the Oj is real: these are the 
conjugate point singularities. When one of the aj is imaginary there are also branch 
points on the imaginary axis. 



7.1 The conformally flat symmetric plane wave 

Consider first the conformally flat symmetric plane wave, with ai = a 2 = a. In this 
case, both polarizations propagate with the same phase velocity and refractive index 
- there is no birefringence. We can therefore set J~{z) = lJ 7 (z). Eq. (|7.4j ) simplifies 
to 

The integral can be evaluated by first rotating the contour to the negative imaginary 
axis and then by direct evaluation, giving a closed-form expression in terms of di- 
gamma functions: 

F+(z) = -zlog^ + + £) - (T ■ (7.6) 
As expected, this is a branched function due to the presence of the logarithm. 

The corresponding function J--{z) defined as 

is given explicitly by 

F.{z) = -z tog(-i) +ziP(l-£)+a . (7.8) 

It satisfies 

T-{z) = -F+(-z) , (7.9) 

17 Notice that this result is slightly different from that quoted in [16,17]. The difference is because 
of the way the overall position of the loop was fixed. In [16, 17] the centre of the loops were fixed 
by hand to be at the origin. In the present work we have not needed to fix the overall position 
of the loops since this is done automatically because the loops are pinned at x(0) to go through 
x, the origin in Rosen coordinates. It turns out that there is a non-trivial Jacobian between these 
prescriptions. 



1 - 



at 



sin at 



(7.5) 



1 - 



at 
sin at 



(7.7) 
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by virtue of the translation invariance of the symmetric plane wave metric, which 
guarantees that the VVM matrix (|7.3|) is a function only of (u — u') and the factor 
[ . . . ] in the integrand of ( |7.5|) and ( [7.7|) is an even function of t. 




Figure 13: The physical sheet for T{z) showing the branch cuts and the simple poles which lie 
on other sheets accessed by moving through the cuts in the direction of the arrows. 

The physical sheet is given by the cut z-plane with the physical T{z) defined as 
the analytic continuation of T+{z) into the lower-half plane and of T~{z) into the 
upper-half plane, that is (see Fig. |T~3D : 




T+{z) = -z\og^ + ztp(l + - a , - 7T < arg z < , 
F_{z) = -zlog(-£) + Z < l f,(l-JL) + <T) < arg z < n . 



(7.10) 



Before considering the discontinuities and the Kramers-Kronig relation, notice 
that in addition to the cuts from the logarithm, T + {z) also has simple poles on 
the negative real axis at z = —2pa, p = 1,2, .. . from the di-gamma functions. 18 
Similarly, T-{z) has poles on the positive real z-axis. These poles are not on the 
physical sheet, as defined above, so do not directly affect the physical refractive index. 
Nevertheless, they encode useful information about the functions J 7 +(z) and F-(z) 
and provide an alternative method of computing the physical discontinuities. 

The full analytic structure of J- + (z) can be understood as follows. First, intro- 
duce a cut-off 6 to regularize the lower limit of the t integral. It is then useful to 

ls ip(x) has simple poles at 0,-1,-2, . . . with residue —1. The di-gamma function also satisfies 
the following identities, to be used later: 

+ z) = lj}{z) + \ , "0(1 — Z ) — 4>{ z ) + 7TCOt7TZ . 
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consider the integral as the sum of two pieces. The first term is 



°°- le dt 



g— iSz 



+ izEi(-i8z) ( 7 - n ) 



= - 5 - Z -\og(-z^)-z( lE -l). 

where the limit 5^0 was taken in the last line. What is interesting is that this 
term accounts for the branched nature of J- + (z); indeed asz^ ze 2m , the exponential 
integral function has a logarithmic branch cut and so the discontinuity of ( |7.11|) is 
—2ttz. The second term, 

poo— it —izt 

Ff>(6,z) = -io 2 \ dt-^— , (7.12) 
J s sm z at 

only has simple poles which can be manifested by expanding the denominator in 
powers of e~ 2tat : 

poo— it °° 

^\5,z) = Aia 2 / dt \2pe- i{z+2pr7)t . (7.13) 
J & P =i 

Performing the t integral 

CO 

^\6, z) = 4a 2 e~ lzS V _^_ e - 2 ^ 
+ v ' ^ z + 2pa 

p=i 

^~ lzS toze-iW , z z 2iaS s (7.14) 

= -^3T - z + 2a 2Fl(1 ' 1 + 5?' 2 + 3?> e ) 



= -- + - log(-^) + z( lE -l)- a + z^(l + ^) . 

Summing the two contributions ( |7.11| ) and (|7.14 ), we see that the divergent terms 
cancel to leave the finite piece ([7.6|). 



Returning to the refractive index, we now perform the integral over ^ in (|67 
and define the physical n(u) on the physical sheet described in Fig. |9| as the analytic 
function found by continuing n + {uo) into the upper-half plane and n-(u) into the 
lower-half plane, that is 

/ \ fn+(u) , < arg^j < rr , 

n{uj) = < ('-15) 
I n_(o>) , — 7r < arg uj < . 

The translation invariance property J~-(z) = — z) implies 

n-(u) = n + (—u>) . (7-16) 
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Also, note that Imn(w) = for real u in this example. With this definition, we also 
see that n(u) is not a real analytic function, i.e. n(u*) ^ n{uj)* . This is because of 
the difference in the functions n + {uj) and n_(u;), and reflects the need for the cuts 
in the z and cu-planes in the definition of T{z) and n(oS). This shows very clearly 
how the geometry, in the form of conjugate points, implies an analytic structure for 
the refractive index which fails to satisfy the usual S'-matrix and dispersion relation 
assumptions. 

Since the symmetric plane wave metric exhibits translation invariance in u, the 
second form of the Kramers-Kronig relation (|6.11| ) should hold in this example. We 
now check this. First, we need the discontinuity Disc T{z) on the positive real z-axis. 
Using standard di-gamma function identities (see footnote 19), we find 19 



Disc T{z) 



7TZ 

nz cot — — inz 
2a 

2mz 



(7.17) 



D iirz/(T 



This can also be found by using the contour in (|6.5|), see Fig. [TT], to evaluate Disc T[z) 
as the sum of the residues of the poles on the real t-axis in the integrand of ( |7.5| ). 
These are double poles at rnr/a, n e Z ^ 0. The fact that the singularities are poles 
rather than branch points is due to the fact that for the conformally flat symmetric 
plane wave conjugate points are simultaneously conjugate for both polarizations. 
The discontinuity in T{z) associated to the series of double poles at t = nn/a is 
then given by (|6.5|) as 



Disc T{z) = — 2ni Res n/K / a 



le 



-izt 



2mz e~ innz/cT 



2mz 



oiirz/a | 



reproducing the result above. 



at \2 



sin at 



(7.18) 



Finally, substituting back into the refractive index formula, we can evaluate the 



19 As a consistency check, we can verify that the sum of the discontinuities on the positive and 
negative axes, viz. 

-2iirz , 



ginz/a ^ g— i-Kz/cr ^ 

which reproduces the discontinuity of the logarithmic cuts in T{z). 
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Kramers-Kronig relation: 



n(0) — n(oc) 



i r°° du 

m J u 

a 



Disc n(u — ie] 



2m 2 
a 



m 2 m 2 j 
a 



rl poo 

/ df[f(l-f)] 2 / dz Disc T(z + ie) 
Jo Jo 



ie 



iir 2 m 2 30 j 
a 1 



dz 



'o 

2mz 



(7.19) 



iw 2 m 2 30 
acr 2 



5/^ _ 1 

2ttz 

oTTz/a ^ 



CO 



n=l 



90nm 2 



where the sum over the residues in the second to last line is evaluated as £(— 1) = — ^ • 
The prescriptions here are crucial because DiscT-^z) has a set of simple poles 
and the integration contour must be defined appropriately. The required definition 
follows from a close examination of (|6.14j ). This shows that "Disc n(u) v is in fact 
n + (uj + ie) — n_{uj — ie) and picks up not just the discontinuity across the cut itself 
but also a contribution from the hidden poles on the unphysical sheet for Re uj > 0. 



We can check this explicitly. For small u, we have from ( p.lOj ) that 



njj(0) = S i:j + 



a 



3607rm 2 



907rm 



i 8 *i 



(7.20) 



For large \u\, the refractive index is determined by expanding ( |7.5| ) for small z. In 
particular, in the limit as \u\ — > oo, we have ( p. 17\ ) 



n(uj) 

plus less singular terms, where 



a 



127TU 



Ho) 



(7.21) 



^(0) 



dt 



5i 



at 



sin at 



00 dt 



'j 



at 



-a . 



sinh at J 

(7.22) 

where we have rotated the contour t — > —it. So we verify that the high-frequency 
limit of the refractive index is indeed 8y as expected, corresponding to a wavefront 
velocities equal to c, and the Kramers-Kronig identity holds (despite the absence of 
a non-vanishing Imn(w)) by virtue of the contribution from Disc n(uj) across the cut 
in the complex a»-plane. 
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The complete frequency dependence of the refractive index can be found by nu- 
merical evaluation of ( [7.5| ) and the result is shown in Fig. |TJ1 Notice that because 
of the sign difference in the one-loop coefficients in ( |3.10| ) and ( |3.9| ) for scalars and 
spinors, the corresponding result for spinor QED is the opposite of this, viz. a super- 
luminal low-frequency phase velocity falling monotonically to c in the high-frequency 
limit. 



7.2 The general symmetric plane wave 



We now extend this analysis to the general symmetric plane wave (|7.4j ). Writing 

Tij{z) = Tj{z)bij, we have 



TAz) 



dt ie 



-izt 



1 



3/2 1/2 

0i °Y 



sin 



3/2 



1 1I 



(7.23) 



and similarly for ^{z) with o\ ^ o~2- With no loss of generality we can assume o\ 
is real, while a 2 can be real or imaginary with <j\ > | cr 2 1 (and we shall assume that 
o\ is positive and argo"2 is either or |). To define a physical sheet, we continue 
J r +j{z) = J-j{z) into the lower-half plane (including the positive real axis) and glue 
it to J z ^j{z) = —J r + j{—z) in the upper-half plane, just as in the conformally flat 
example. 



To demonstrate the new physics arising here, we numerically calculate the re- 
fractive index for the Ricci flat symmetric plane wave metric, G\ — a and cr 2 = iu. 
This displays gravitational birefringence, in that the two polarizations move with 
different refractive indices. Moreover, in this case, the refractive index also develops 
an imaginary part, which would not be seen in the low-frequency expansion based 
on the effective action. The real and imaginary parts of the refractive indices are 
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plotted in Fig. [15]. Remarkably, and unlike the usual case of dispersive scattering 
where Im n{u) > 0, we find that here the imaginary part of the refractive index is 
negative. 
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(b) 

Figure 15: (a) Reni(w) — 1 and (b) Imtiifw) (i = 1 continuous, i — 2 dashed) in units of 
aa 2 1 (irm 2 ) , as a function of log ua/m 2 for the Ricci flat symmetric plane wave. 



Although we do not have a complete expression for the refractive index in this 
case in terms of elementary functions, we can still get a very accurate approximation 
using analytic techniques. First of all, we rotate the contour by taking t — > —it: 



T{z) 



dt e 



-zt 
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sin s at sinh 2 s at 



(7.24) 



where s = \ or | depending on the polarization. The integration contour lies on top 
of the branch points at t = nn/a, n — 1, 2, . . .. Notice that the integrand is real for 
0<t<7i,27i<t< 3tt, etc., and imaginary for 7r < t < 2n, 3n < t < Att, etc. 
Since the integrand is falling off exponentially like e -( z +( 2 - s ) a ) t we can approximate 
the imaginary part by expanding around the first branch point t = n/a; first, taking 
s 
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2 3/ 2(T 3/2 
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(7.25) 
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and then for s 
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2' 
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2(2e-")5vW /2 e 



l/2.-«7r/<7. 



Z 



a/2 



20 The integral here appears to be singular at the lower limit. However, in reality the contour 
jumps over the branch point and this regularizes the integral in a way which is equivalent to taking 



J oo ^e-7i 3/2 -r(-l/2) = -2 v ^. 
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From these expressions, we can get a further approximation of the refractive 
index itself valid for low frequency, by evaluating the £ integral around the saddle- 
point of the exponential factor e~ Z7r /°" = exp — Tim 2 / {2uja^(l — £)) which occurs at 
£ = ~. This gives the leading low frequency behaviour as 

2 

ImniM = — ^— (2e-") 3/2 v^Fe- W/(w<T) + ■ • ■ (7.27) 
327im 2 

and 

\mn 2 {uo) = - — (2e- n ) 1/2 ^e- 2 ™ 2/(uJa) + • • • , (7.28) 

87TCo> 

which accurately reproduces the numerical evaluation in Fig. 

To try and understand the origin of this unusual dispersive behaviour, we can 
follow the same logic as for the conformally flat example to deduce the analytic 
structure of !F + j (z) , since ultimately the sign of Imn(w) is determined by the location 
of branch points on the unphysical sheet of n{oj). The idea is to introduce a cut-off 
5 on the lower limit of the integral and consider the contribution from the two terms 
in the integrand. The contribution (|7.11|) remains the same, whereas the second 
contribution is now 



roo-ie 3/2 1/2 

J?\6, z) = - dtie~" , 3/2 \ 2 1/2 I • ( 7 -29) 

Jo sin 6 ' ait sin 1 a 2 t 

We can expand the integrand in terms of e~ 2lr7lt and e -2 "^*, which is a convergent 
expansion along the integration contour. (Notice that this is the expansion which is 
consistent with our choice of argo"2 to be either or |.) Performing the t integral 
on the terms in the double expansion gives 



3 1 
3\ /1\ e -i(2m+^)o-i5-i(2n+7j)o-J5 



While we cannot sum this in closed form, we know that apart from a singular term 
which cancels that in J z + j{z) ( 7.11 ), it is a holomorphic function with simple poles 



at 

z = -(2m+§)<7i-(2n + i)<72 , m,n eZ>0. (7.31) 

In particular, in the case when a 2 is imaginary there are poles in the upper-half 
plane. A similar story holds for the other polarization state. 

To conclude, J z + j{z) has a branch point at z = coming from the zlogz term 
in ( [7.11 ) along with a set of simple poles which lie in the region Rez < 0, Imz > 0. 



In particular since they lie in the upper-half of the z plane they give rise to branch 
points in the lower-half of the u plane (on an unphysical sheet) and are therefore in 
the causally safe region. 
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The situation for U2(uj), where the polarization lies in the direction of diverging 
geodesies in the null congruence, is therefore quite similar to the simple single-pole 
refractive index model discussed in Section 6.4 in that we find a low- frequency su- 
perluminal phase velocity together with branch points in the negative imaginary half 
cu-plane. The resulting Imw2(w) < implies an amplification of photon propagation 
through this background spacetime, centred on a characteristic frequency of order 
^R/m 4 ~ 1. In contrast, ni(u>), where the polarization lies in the direction of con- 
verging geodesies, shows similar behaviour to the conformally flat example, with a 
subluminal phase velocity and only a very small, though still negative, imaginary 
part Imn 1 (co>) < 0. 



8. Example 2: weak gravitational wave 



We now consider an example of a time-dependent background, the weak gravitational 
wave, which displays gravitational birefringence and dispersion with both positive 
and negative imaginary parts for the refractive index. 

The spacetime metric for a weak gravitational wave takes the following form in 
Rosen coordinates: 

ds 2 = -2dudV+ (1 + e cos vu)dY l dY l + (l - e cos vu)dY 2 dY 2 . (8.1) 

Here, and in the following, e is small and we work to linear order. The transformation 
to Brinkmann coordinates is achieved via the zweibein 

Ej(u)=S i a (l-{-iylcoBvu) , (8.2) 

to give 

■ ev 2 

ds 2 = —2du dv + (— l)- 7 -^- cos vu y^y^ du 2 + dyi dy^ . (8.3) 

The equation for the Jacobi fields is 

d 2 y^(u) ev 2 

d J } = T— cos vuy>(u) , (8.4) 

with ± corresponding to j = 1 and j = 2, respectively. This can easily be solved 
perturbatively in e, with solution to linear order 

£ 6C2 

y (±) (u) — c\ + uc 2 ± -(ci + UC2) cos vu =F — sin vu . (8.5) 
2 v 

Solving the Jacobi equation ( |5.6j ) with the boundary condition (|5.7|), we now find 

the eigenvalues of A 

A {±) {u,u') =u-u' 

,£(«-«')/ A ^ I • a (8 ' 6) 

± y cos vu + cos vu) =F — [ sin vu — sin vu ) , 
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which determines the eigenvalues of the Van-Vleck Morette matrix: 



A (±) (u,u') = 1 =F - (cos vu + cos vu) ± — 

2 v[u — u') 



sin vu — sm vu 



(8.7) 



The refractive index is given by fl4.42|) with J-[u] z) = diagjF (±) (-u; z), where 



F^Hmz) = ±e 



dt 

T 2 



ie 



-izt 



-(cos vu + cos v{u — t)) ^ ( s * n vu ~ s * n u ( u ~ ^)) 



=l=e \fi(z) cos vu + f2{z) sin vu] 



where 



fi(z) = — [2v(l + log z) + (z - v) \og(z -u) — (z + v) \og(z + V 
f2(z) = — [v 2 + 2z 2 log z + z{v - z) \og{z - v) - z{z + v) \og{z + v)\ 



■9) 



The functions fj(z) have branch points at 0, oo and z = ±v and this means 
that rij(u;uj) will have branch points at ±oo and ±2m 2 /v. In particular, the 
branch points at ±2m 2 / v are points of non-analyticity of the refractive index. This 
non-analyticity manifests itself by the fact that lmfi(z) and Re/2(^) are zero for 
z G R > v, while for z G R < v, 



lmf 1 (z)=Ref 2 (z) 



7TZ{V — Z) 

Av 



.10) 



It is then a simple matter to extract the low frequency expansion of the refractive 
index: 

aev 2 r 1 1 /'ujv\ 2 



n (±, (u:uj) = 1 ± 



m 2 ir 
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.aev 
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while at high frequencies, 



nj{u;u) 



- , - aei/ ■ 

1 ± i sm vu + 



cos vu 



sm z/m 



(8.11) 



(8.12) 



The full form of the frequency dependence of the real and imaginary parts of the 



refractive index is plotted numerically in Fig. [16], evaluated at a fixed point on the 
photon trajectory. For the first polarization, this shows a conventional dispersion 
for Ren(w), with a single characteristic frequency uj 2 R/m i ~ 1, together with the 
corresponding Imn(w) > 0. Once more, however, the second polarization is superlu- 
minal at low frequencies and has Imn(w) < 0, indicating amplification rather than 
dispersive scattering. The roles of the two polarizations of course change along the 
photon trajectory through the background gravitational wave. 
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Re ii- 1 Im n 




Figure 16: (a) Ren (±) — 1 and (b) Imn (±) for u = 0.2, m = v = 1 plotted as a function of logo; 
in units of ea for both polarizations. The point of non-analyticity at uj = 2 is quite clear. (Note 
that the fact that the polarizations do not quite give mirror images is an artifact of the numerical 
approximation. 



9. Green Functions 



So far, we have been considering the one-loop correction to photon modes that come 
in from past infinity and propagate out to future infinity. As we have already pointed 
out, this is, strictly-speaking, only consistent if the space becomes flat in those limits, 
otherwise the one-loop correction to the mode becomes large undermining pertur- 
bation theory. A local way to investigate causality involves specifying some initial 
data on a Cauchy surface and seeing whether it propagates causally. This avoids the 
problem of having modes come in from the infinite past. Such initial value problems 
lead to an investigation of the Green functions. 

For a general spacetime, it is not possible to construct the complete Green func- 
tions due to the fact that we can only construct the modes in the eikonal limit 
u ^> VR. However, if we are interested in the one-loop correction to a Green 
function G(x,x'), in the neighbourhood of the component of the light cone with 
V — V — and Y a = Y' a = — and this will teach us about the one-loop correc- 
tion to the causal structure — then it is consistent to replace the full metric by the 
Penrose limit of the null geodesic which goes through (V = 0, Y a = 0). In this way, 
we need not work in the eikonal approximations since the modes ( |2.2| ) are exact in a 
plane-wave spacetime. Once we have taken the Penrose limit, then it is possible to 
calculate the one- loop correction to the Green functions exactly. (See Fig.|l7D. 

In order to construct the Green functions we need a complete set of on-shell 
modes. The most immediate problem is that the general plane wave spacetime does 
not admit a set of Cauchy surfaces. However, we will follow [39] and use the null 
surfaces u = const, to define the canonical structure, a choice which is sufficient for 
our purposes. We now search for the complete set of on-shell modes with respect to 
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u 



V 



Figure 17: The one- loop correction to the Green functions can be calculated for two points 
x' = (u',0,0) and x = (u,V,Y a ) for small V and Y a , so in the neighbourhood of the light cone, 
since the full metric may be approximated by the Penrose limit as indicated. 



the inner-product defined on the ersatz Cauchy surfaces: 

(A, A') = ~ J dVd 2 Y Vdet^(w) [A ■ d v A'^ - A /f • d v A) . (9.1) 

The modes with 

A (i)u (x) = A m {x) = , A (l)a (x) = A{u)E ia {u)e-^ v , (9.2) 

that we constructed in Section 3 are clearly on-shell, but are not the most general 
set of modes. A complete set of gauge-fixed on-shell modes can be constructed by 
taking a more general eikonal phase and polarization: 

A (i)lM (x) = e (i)M (x) exp -i [uV + p a Y a + V(u)] , (9.3) 

where the eikonal equation (|2.4j ) determines 

dlp(u) 1 v-,_ Xl ,iab 

-± l = ^[C '(«)] Ml, (9-4) 

implying, for later use, 



tf(ti) -#*0 = ^^[A-VX)] a W& • (9.5) 
The gauge- fixed polarization vectors now pick up an additional component: 

A J7 ( \ A PaEi a (u) 

00 

while the scalar amplitude remains as in (|4.31 ). Notice, in Brinkmann coordinates 
the polarization vector is particularly simple with i^j = S^j and e^) u as above. The 
modes are split into the positive/negative frequency on-shell modes as according to 
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whether uo ^ and we will define the 3-momentum p = (ou,p a ). One easily finds that 
the inner-product on the modes is 



{A (l) (p),A u 0)) = 2(27r)W 8 >(jf- 



(9.7) 



All the various propagators can be constructed from these modes. We begin by 
constructing the Wightman functions 21 



(9.1 



where, as indicated, the integral over uj extends over to ±oo for G^, respectively. 
The p a integrals are Gaussian and hence easily performed: 
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(9.9) 



where vo^x^x') has the following non-zero components in Brinkmann coordinates: 



We therefore find 
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iu<r(x, x n 



u — u' 



with 



a(x,x') = (u - u')(V - V) + \A ab {u, u'){Y - Y') a (Y - Y') 



(9.12) 



The Feynman propagator iG F ^ lI/ (x,x') = (0\T(A ll (x)A u (x')^ |0), is given by 

iG F ^(x, x') = 6{u - u')G+ v {x, x') + 9{u' - u)G~ u (x, x') . (9.13) 



21 



We use the notation of Birrel and Davies [40] . 
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However, in the present context we are more interested in the causal Green functions. 
In particular, the Pauli- Jordan, or Schwinger, function 



iG^(x,x') = (0|[A M (x),A,Or')]|0) 



(9.14) 



is 



iG^v^x, x ) G (x 1 x ) G (x, x ) 



1 w^ u {x,x' 



,A /-00 



2(2tt) 



u — w 



du exp 
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(9.15) 



— G7„^(x, rr')\/det A(w, m')5(<t(x, a;')) 
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From this the retarded and advanced Green functions may be extracted via 
G Rflv {x,x') = -0{u- u^G^x^x') , G A/JiV (x,x') = 9(u' - u)G^(x, x') , (9.16) 

The causal properties of the Pauli- Jordan function are manifest in the last line 
of ( 19.151) . Gn U (x, x') has support only if x lies on the forward or backward light cone 
of x'. In particular, for Gn fMll (x,x / ) the support is on the forward light cone and for 
Ga^u(%, x') it is on the backward light cone. This is precisely what is to be expected 
for the causality properties of the Green functions of massless quanta. 



The one-loop correction to the Feynman propagator is given by the usual ex- 
pression 



Gf7( x > x ') = ~ j d ' xd ' x ' {g(u)g(u')) 1/2 G F ^(x,x)U 



l-loop^ ~/ 
crp 



x, x )G 



fj,{x i X 



(9.17) 

As explained above, the causal properties of a theory are not manifested directly in 
the Feynman propagator, which receives contributions both inside and outside the 
light cone. In order to address the causal structure we need to calculate the one-loop 
correction to the Pauli- Jordan function, or retarded and advanced Green functions. 
However, these can extracted from (|9.17|) using ( |9.13| ) and ( |9.15|) . 



The key result is a generalization of the calculation of Section 4 with the more 
general on-shell modes Q9.3Q: 
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where the ellipsis indicates additional terms that do not depend on the curvature, 
i.e. are needed to have the correct flat space limit. 



The strategy for calculating (|9.17 ) is to write the tree-level Feynman propagators 
in terms of G using ( |9.13| ). Then we write G^ 1 in terms of the on-shell modes, 
as in ( |9.8| ). Once this has been done we can use ( |9 . 1 S[ ) . The key point is that 
( p,18| ) conserves the "momentum" p and so the contributions schematically of the 
form G + UG~ and G~UG + vanish, leaving the non-vanishing contributions G ± HG ± 
which are immediately identified as the one- loop corrections to G^. Notice that the 
step functions that are present in ( |9.13|) mean that u and u' are constrained to be 
u > u > u' > u' and u < u < u' < u', respectively, for G 1 ' 100 ^. It is then convenient 
to change variables from {u, u') to (u, t), where t — u — u' . Putting all this together, 
we have 



2ia 
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x Aij(u, u — t)A/detA(tt, u — t) 

A {i)fl (p;x)A {j)u (p;x'y + 



(9.19) 
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which are initially valid for u ^ u', respectively, but which can be extended to all u 
and vl by analytic continuation. 



The p a integrals in ( |9.19|) are identical to (|9.9| ) and so the former becomes 
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(9.20) 
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where we have just displayed the components with spacetime indices in the two- 
dimensional polarization subspace. 



As happened at tree level, the Pauli- Jordan function ( |9. 15| ) is given by ( |9.20|) by 
extending the u integral from — oo to +oo: 



£~<l-loop^ 
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(9.21) 
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Notice that in the limit v! — > —oo with fixed it, the expression can be written in 
terms of the refractive index. In particular, for x = (u, V, 0, 0) and x' = (u', 0, 0, 0) 
in the limit v! — > — oo, the Pauli- Jordan function (or, since u > u', the retarded 
propagator) is 

/u /*oo 
du / dujn i3 (u;uj)e- iujV . (9.22) 
-oo J — oo 

Since n(u; u>) in analytic in the upper-half plane, the retarded propagator vanishes 
when V < 0, i.e. outside the light cone. This confirms that even in the presence 
of the novel dispersion relations and superluminal phase velocities described here, 
causality is maintained with advanced, retarded and Pauli- Jordan Green functions 
displaying the necessary light-cone support. 



A remarkable feature of ( |9.21D is that we can rewrite it in a manifestly causal 
form in terms of the Pauli- Jordan function G(m 2 ;x,x') of a massive scalar particle 
(see eq.([OD): 

^ lo ° P 0r,*') = - fdu r U % I' dit{\~t) 

71 Ju ' Jo /° 9 v (9.23) 

/ 777, t \ 

x A^u - t)^detA(u,u - t) G — — jr]x,x' +■•• . 

As in previous formulae, the ellipsis represent terms that do not depend on the 
curvature. This last expression makes the causal structure completely manifest. In 
particular, G(m 2 ;x,x f ) has support only inside, or on, the light cone and so at the 
one-loop level the commutator of two photon fields receives contributions from inside 
the light cone. However, causality is maintained because the one-loop correction 
Gl~ loop (x, x') still vanishes outside the light cone. 



10. Summary and Conclusions 



In this paper, we have analyzed the effect of vacuum polarization in QED on the 
propagation of photons through a curved spacetime background. This problem is 
of potentially fundamental significance because of the discovery that quantum loop 
effects can induce a superluminal phase velocity, raising the question of how, or 
whether, this can be reconciled with causality. We have resolved this issue through 
an explicit computation of the full frequency dependence of the refractive index for 
QED in curved spacetime, showing that the wavefront velocity, which is the speed 
of light relevant for causality, is indeed c. This is, however, only possible because 
a number of generally assumed properties of QFT and 5*-matrix theory, including 
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the familiar form of the Kramers- Kronig dispersion relation, do not hold in curved 
spacetime due to a novel analytic structure of the Green functions and scattering 
amplitudes. 

The key insight which makes this analysis possible for general spacetimes is the 
realization, inspired by the worldline formalism of QFT, that to leading order in 
the curvature, the quantum contributions to photon propagation are determined by 
the geometry of geodesic deviation, i.e. by fluctuations around the null geodesic 
describing the photon's classical trajectory. This geometry is encoded in the Penrose 
limit of the original spacetime, simplifying the problem of photon propagation in 
general backgrounds to that of their plane wave limits. 

This allowed us to derive a compact expression for the complete refractive in- 
dex n{uj) of a curved spacetime entirely in terms of the Van Vleck-Morette matrix 
for its Penrose limit. In this form, the novel analytic structure we have discovered 
becomes manifest. It is related to the occurrence of singularities in the VVM matrix 
corresponding to the existence of conjugate points in the null geodesic congruence de- 
scribing photon propagation. The existence of conjugate points is a generic property 
of geodesic congruences, related to the Raychoudhuri equations and enforced by the 
null energy condition. This geometrical origin shows that the type of unconventional 
analytic structure which we find here for n(u), notably the loss of real analyticity, 
will also occur in more general scattering amplitudes in QFT in curved spacetime. 

The analytic structure of the refractive index and its implications for causality 
were described in detail in Section 6. The VVM singularities mean it is necessary to 
define n(uj) on a physical sheet, pasting together analytic functions on opposite sides 
of cuts running along the real w-axis. This construction banishes dangerous branch- 
point singularities to the unphysical sheets, leaving the physical refractive index 
with the crucial property of analyticity in the upper-half complex u;-plane required 
for causality. We also investigated causality at the level of the Green functions 
themselves, demonstrating explicitly how the retarded, advanced and commutator 
functions have support on the relevant parts of the fundamental light cone. 

The general theory was illustrated in two examples - the symmetric plane wave 
and the weak gravitational wave background. As well as demonstrating explicitly 
how the analytic structure of the refractive index arises and how the Kramers-Kronig 
dispersion relation is realized in a way compatible with causality, these examples 
introduced a further surprise, viz. the occurrence of negative imaginary parts for the 
refractive index. Such behaviour is generally associated with the amplification of 
a light wave as it passes through an optical medium exhibiting gain, rather than 
the usual dispersive scattering associated with the production of e + e~ pairs. It 
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appears that in this case quantum loop effects in the curved spacetime background 
are responsible for the emission of photons. 



In addition to clarifying the mechanism responsible for inducing Imn(w) < and 
understanding the relation to the optical theorem, there are many other directions in 
which this work can be extended. An immediate task is to generalize the fundamental 
result ([4.41 ) for the refractive index to spinor QED and other QFTs. This requires 



some further technical developments, but the discussion of analyticity is unlikely to 
change significantly. It is also interesting to extend these results to cosmological and 
black hole spacetimes and study the effect on photon propagation of horizons and 
singularities. The phase velocity derived from the low-energy effective action already 
displays special simplifications at a horizon [5,41], while the Penrose limit of black 
hole spacetimes near the singularity exhibits an interesting universality [42] that will 
be inherited by the refractive index. 

Finally, it is important to extend this analysis from photon propagation and 
two-point Green functions to higher-point scattering amplitudes. In view of the geo- 
metrical origin of the novel analyticity structure discovered here, it seems inevitable 
that many of the conventional analyticity properties of scattering amplitudes, which 
underlie all the usual relations of 5-matrix theory, will also be radically changed 
in curved spacetime. In particular, a study of the analytic structure of Planck en- 
ergy scattering may well have far-reaching implications not only for QFT in curved 
spacetime but also for string theory and perhaps even quantum gravity itself. 
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